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^ ■ I. INTRODUCTION 

> ■ 

The study of strongly correlated electrons continues to receive a lot of interest due to applications in condensed 
matter physics. Some of the well known models are the Hubbard and t-J models [|l| and generalizations such as the 
EKS model @. 

f*^ Another example of a model describing strongly correlated electrons is the the supersymmetric U model. This 

0^ model was first introduced in |^ and was shown to be integrablc via the Quantum Inverse Scattering Method (QISM) 
0^ by demonstrating that the model could be obtained from an _R-matrix which is invariant with respect to the Lie 

superalgebra gl(2\l). The Bethe-ansatz equations for the model were obtained in ||^-^. Subsequently, an anisotropic 
generalization was presented in |9| which was also shown to be integrablc through use of an i?-matrix derived from a 
representation of the quantum superalgebra Uq{gl{2\l)) pO| . 

■ The anisotropic supersymmetric U model describes a system of correlated electrons and generalizes the Hubbard 
Q . model in the sense that as well as the presence of the Hubbard on site (Coulomb) interaction there are additional 
^ ' correlated hopping and pair hopping terms. The model acts on the unrestricted 4'^-dimensional electronic Hilbert 

space (^n=i^^ where k is the lattice length. This means that double occupancy of sites is allowed and distinguishes 
this model from the anisotropic t — J model which shares the same supersymmetry algebra Uq{gl{2\l)). The model 
^ . , contains one free parameter U, (the Hubbard interaction parameter) which arises from the one-parameter family of 

■ inequivalent typical four-dimensional irreducible representations of the Uq[gl{2\l)], and another which arises from the 
deformation parameter q. 

Bethe ansatz solutions for the anisotropic model with periodic boundary conditions have been studied P, p^ , p^ , 
however for this case there is no quantum superalgebra symmetry. In p^-|l8|| some quantum algebra invariant 
integrablc closed chains derived from an i?-matrix associated with the Hecke algebra were introduced and investigated. 
It was subsequently shown ||l^] that a general prescription for the construction of integrablc systems with periodic 
boundary conditions and quantum algebra invariance existed which could then be applied to higher spin models such 
as the spin 1 XXZ Heisenberg chain po[ |. 

In the present article we further develop this method by considering the graded case to derive the Hamiltonian of 
the anisotropic supersymmetric U model with quantum supersymmetry on the closed chain. We will adopt a nested 
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algebraic Bethe ansatz to solve the model and this will be presented in detail in Section 3. Also the energy of the 
Hamiltonian will be given. 

II. QUANTUM ALGEBRA INVARIANT HAMILTONIAN FOR THE SUPERSYMMETRIC U MODEL 

The following notation will be adopted. Electrons on a lattice are described by canonical Fermi operators Ci^a and 
c| ^ satisfying the anti-commutation relations given by {cl ^, Cj^r} = SijSar, where i,j = 1,2,.., A: and tr, r =t, J. • The 
operator Ci,a (c|o-) annihilates (creates) an electron of spin a at site i, which implies that the Fock vacuum |0 > 
satisfies Ci^a\0 >— 0. At a given lattice site i there are four possible electronic states: 

io>, it>4=4tIo>' \^>^^clJo>, I ti>.= 4j^4^|o> . 

By rii^cr = cl ^Ci^cr we denote the number operator for electrons with spin a on site i, and we write = + ni^^. 
The local Hamiltonian for this model is ||^ 

Hi(i+i) = ~^{cl^Ci+ia + h.c.)exp 

a 

+ ^Uni-fnii + Uni+i-fni+ii + U{cl^cl^Ci+iiCi+i-^ +h.c.) , (1) 
where i labels the sites and 

[/ = e[2e~^cosh - cosh 7)]^, e = ±1. 

This Hamiltonian may be obtained from the i?-matrix of a one-parameter family of four-dimensional representations 
of Uq[gl{2\l)], which is afforded by the module W with highest weight (0, Ola). The details of this construction may 
be found in pO[ |. 

The Hamiltonian (|^) may be modified to ensure quantum superalgebra invariance by adapting the general con- 
struction presented in ||l9| . We can write 

4=1 

where the boundary term is given by 
with 

G = R2iR^i---R^i, k the lattice length. 

Above, i?^ is the constant i?-matrix obtained as the zero spectral parameter limit of the Yang-Baxter equation solution 
associated with the model pO| . These operators act in the quantum space and the closed boundary conditions of the 
model may be explained by the relations 

GHi^i+i ^ Hi+is+2G, i — 1,2, ...,k — 2, GHik — H12G. 

The quantum supersymmetry of the Hamiltonian is a result of the intertwining properties of the matrices R. 

III. NESTED ALGEBRAIC BETHE ANSATZ 

We present the nested algebraic Bethe ansatz for the above Hamiltonian by extending the methods presented in 
p9| , p0| to treat quantum group invariant closed higher spin chains to the graded case. We begin with the _R-matrix 
satisfying the Yang-Baxter equation constructed directly from a solution of the twisted representation as given in . 

The Yang-Baxter Equation may be written as the operator equation: 
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acting on the spaces V ®V where V is the vector module and W is the four-dimcnsional module associated with 
the one-parameter family of minimal typical representations. Greek indices are used to label the matrix spaces, that 
is the first two spaces and Roman indices label the quantum space, which is the third space. The quantum space 
represents the Hilbert space over a site on the one-dimensional lattice. The „i,i?-matrix acts in the matrix space and 
it is between the two matrix spaces that the graded tensor product acts. 
The ^t,i?-matrix acts on F (g) F and has the form plj, 
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where A, E and C depend on the spectral parameter as follows Aix) = ^i^] 
t,„i?-matrix satisfies the Yang-Baxter equation 

Ri2ix/y)Ri3ix)R23{y) = R2MRMRi2{x/y). 

By construction, these i?-matrices also satisfy the generalized Cherednik reficction property 

Rl?P'{^)R~%'\^ly) = Rl?,,{y)R'%''{ilxl 



and C{x\ 



The 



and crossing unitarity pA 



R'K^0,{xOK'4R-'r€f{x)K-'^' = 
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where sti denotes matrix supertransposition in the first space and K is the crossing matrix given below. 
It will be necessary to rewrite the matrix in terms of constant matrices yyR^ and ^vR^ that is. 
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where yyR~^ {wR ) corresponds to the leading term in the limit as x — > 00 {x —^ 0). 

The ^tui?-matrix was constructed in [Q in the V ®W space and has the following form 
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where the dependence of these elements on the spectral parameter is given by 



J{x) 



{xq-"- 



-a-2 



Y{x) = J{x){D + B) + 
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L{x) = 



^^([a + l]J(x) + l), M{x)^F^DJ{x) + j^, N{x) = ^^{J{x) + [a + 1]), 
[a + 2\ [a + 2\ [a + 2\ 



Q{x) = (qB - Dq-^)J{x) 



2] 



, Q'{x) - [q-^B ~ qD)J{x) 



2]' 



S{x) 
T(x) 



[a + 2] 



^-(a+3)/2 _ q('^+i)/2FDJ{x), S'ix) 





-2]' 




+ 1] 



[a + 2] 
P{x) = q(''+^y^FDJ{x) 



2] 



v-^; - [a + 2] 



with constants 



D = 



a + 2]iq + q-^)-' 



F 



, B = l/{q + q and [a] 



+ 2]^ 



q" - g-" 
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The i?-niatrix as weh as satisfying the Yang Baxter relation (|^), also satisfies the generalized Cherednick reflection 
property and crossing unitarity (^). 

We now introduce an auxiliary doubled monodromy matrix 
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acting on y (g) W®^ . Above vwR+ represents the leading term in the matrix y^R^^{x) for the limit as a; — + oo. 
Represent the doubled monodromy matrix in the following way: 



03 



.UZ{x) 



vwU-i [x) vwU2 [x) vwU2 (x) 



(5) 



u:Ui{x^ vwU2{x^ vwU^i^X^ 

It may be shown that this monodromy matrix satisfies the following Yang-Baxter relation 

vvRoilfiiliy/^) vwU^^i^{x) yyR^^^^g^ vwUg^^{y) = vwU^^i^ {y) ^yR^g^^^ i,tuC/;^2^(x) (6) 
depicted graphically below. 
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The occurance of the constant matrices vvR-\- will greatly simplify the calculations of the algebraic Bethe ansatz. 

In order to perform the nested algebraic Bethe ansatz (NABA) we define an auxiliary transfer matrix as the (super) 
Markov trace of the monodromy matrix, that is, 
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vwTj-x (y) 



1 J2 




(7) 



where 




(8) 



Therefore the vwT{y) form a one-parameter family of commuting operators and it may be shown that they commute 
with the transfer matrix, wwT{y) of the Hamihonian (|l|). This means that they have a common set of eigenvectors. 
We find that 
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Take the lowest weight state as a reference state (pseudo-vacuum) in W, which we denote as |0 >i. Then |0 >= 
(8)^=1 10 >i and we find the action of the doubled monodromy matrix on this reference state to be given by 



,U{x)k\0>-~ 



I{xf 
I{xf I |0 >, 



(10) 



where 
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We construct a set of eigenstates of the transfer matrix using the technique of the NABA. The creation operators 
are vwUi{x), «toJ7|(a;) due to the choice of reference state. Thus we use the following for the ansatz for the eigenstates 
of vwT{y): 



Vl/ = „^C/,\(xi) ,^Ul_{x2)....^Ul{Xr) *Si\|0>, 

where indices have values 1 or 2. We seek a solution of the eigenvalue equation 

vwTiy)"^ = ™A(y)*. 

The action of these states is determined by the monodromy matrix and the relations 
the construction of the NABA are 



(11) 
(12) 

The relations necessary for 
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with the indices taking values of 1 and 2. It can be seen that this i?-matrix r{y) also fulfills a Yang-Baxter equation 
and can be identified with the i?-matrix of the quantum spin ^ Heisenberg (XXZ) model. The action of the ansatz 
( prj) on the diagonal elements of the monodromy matrix ^ is given by 



q- fJ^Eiy/x,) 



+ u.t. 
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j=i ^^jlV) 

where 

.»T«(jy)-g-i ,^(7«l(y) + g „^f/W^(y). (15) 

In order that the eigenvalue problem ( |l^ ) is satisfied, it is necessary to solve a new eigenvalue problem for the 
nesting as follows: 

.«,TW(y)*W=AW(y,{y,})vI/W, 

where 

^^'^ = v^U^'^liyi) v^U^'^l{y2)-v^U^'^l{yrn)\0 >(i) . (16) 
The second level reference state is given by |0 >^^-'= ®[^i|2 >i. We represent the nested monodromy matrix as 

The action of the nested monodromy matrix vwUm\y) on the second level reference state is 

r 

„^t/«}(y)|0>W =q^l[E{x,/y)\0>^'\ 

r 

„„t/(i)2(y)|0>(i) ^l[A{x,/y)\0>^''> . 

(18) 

The action of vwT^^Hu) '^^ the ansatz ( |l6| ) is computed similar to the first level case from the relations (||). We 
obtain 



.C/(i)^(,)vI/(i) =q-f[ jQ A(.,/y)v,(i) + u.t. 



The eigenvalues for the auxilliary transfer matrix, vwT{y) are found to be 



t. 



(19) 



..A(.) - .-nrf- n - n Hi ft 1^ - (-1) v- ft (») 

provided the "unwanted terms" cancel. The cancellation of these terms lead to the Bethe ansatz equations obtained 
by eliminating the poles of the eigenvalues ( pO| ) 

in 

g'»/(x,)^-^n f / = l,...,r. (21) 

q Vj- xiq 

Associated with these solutions, the energies of the Hamiltonian are given by 

Y (q"/2a;T'/'-g-"/2xy')(g-"/2-ix7'/'-g"/2+ia;'/') 
This expression reduces to the normal periodic case |5| in the rational limit as q — > 1. 
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IV. CONCLUSION 



In this work, we have constructed a quantum algebra invariant supersymmetric U model on a closed lattice and 
derived the Bethe ansatz equations. Notice that in the Bethe ansatz equations (pT| ) the presence of "g" terms in 
comparison with the correspond ing equations for the usual periodic boundary conditions [ p^ . In fact, this feature 
also appeared in other models jr^ , p^ i"8|j25| ] and seems to be a peculiarity of quantum-group-invariant closed spin 
chains. In the limit for g — > 1, the usual Bethe ansatz equations for the periodic chain |^] are recovered. 

An appealing direction for further study of the present closed supersymmetric U model with C/q[(7/(2|l)] invariance 
would be to investigate its thermodynamic properties. In particular the partition function in the finite size scaling 
limit which can be used to derive the operator content of the related statistical models. 
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